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Abstract 



Using the Newtonian approximation, we study rotating compact 
bosonic objects. The equations which describe stationary states with 
non-zero angular momentum are constructed and some numerical re- 
sults are presented as examples. Limits on the applicability of the 
Newtonian approximation are discussed. 
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1 Introduction 



Boson stars are well known gravitionally bound states of complex scalar 
fields. These objects were first studied by Ruffini and Bonazolla [1], and, 
since then, a large number of papers on this subject have been published 
[2], including three recent reviews [3]. As pointed out by Ferrell and Gleiser 
[4], we have two distinct motivations to study boson stars. On one hand, 
they make a good laboratory to study compact self gravitating objects and 
to explore the differences and similarities between then and the usual stars. 
On the other hand, considering that scalar fields play an important role in 
theories of fundamental forces, it is reasonable to address the question of 
formation and stability of bosonic compact objects, including also the study 
of mixed boson-fermion objects [5], which may be even more probable to find 
in nature. 

Concerning the analysis of stable solutions of boson stars, most of the 
work done is restricted to spherically symmetric configurations, including the 
ground state and excitations which are also spherically symmetric. Among 
the exceptions, there is the work of Ferrell and Gleiser [4], which analyses 
the emission of gravitational radiation by boson stars in the Newtonian ap- 
proximation. The compact object is supposed to decay from an excited state 
with non zero angular momentum to the ground state with the emission of 
gravitational waves. The excited state is treated as a pertubation in which 
the excited bosonic particles, carrying the angular momentum, move in the 
background potential of the spherically symmetric state. 

Recently, Kobayashi, Kasai and Futamase [6] consider the slow rotation 
of a relativistic bosonic star. They look for slowly rotating solutions which 
are similiar to the rotation of conventional objects, following the approach of 
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Hartle [7]. They conclude with a negative result, namely that the relativistic 
boson star has no stationary solutions with slow rotation. However, the 
possibility of rapid rotations, which can not be treated pertubatively, is not 
excluded. 

The failure to describe rotational states pertubatively is itself instructive. 
First order pertubation theory is well suited to describe pertubative states 
which can be obtained by continuos deformations of the ground state and, 
to keep it first order, we must also require that the pertubation parameter is 
"small" in some sense. If the excited rotational states we are searching for 
form a discrete set, they can not be obtained by pertubation theory, even if 
they do not rotate at relativistic speeds. And up to now, we do not known 
whether this compact objects always rotate at relativistic speeds or not. 

In this paper, we return to the problem of rotating boson stars and, in 
an attempt to avoid the pertubative approach, we consider the Newtonian 
approximation, which is known to be valid for boson stars provided that the 
central density or the total mass of the star is not higher than a certain critical 
limit [4]. We assume the hypothesis that, at least for the first excited states, 
the rotational effect may be well described by the non-relativistic theory. 
This is analogous to the excited states of atoms which are not obtained as 
pertubations of the ground state but, even so, are well described by the 
non-relativistic theory. Under this assumption, the boson star is analysed, 
and we search for stationary solutions with non-zero angular momentum, 
following the approach of [4], with one major difference: the excited sates of 
the scalar field and their gravitational potential are computed simultaneously 
in a coupled system of equations, without any reference to a fixed spherically 
symmetric background. The resulting solutions describe stationary rotational 
states which are not pertubative deformations of the ground state. At the 
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end, we check the validity of the Newtonian approximation and we identify 
the limits which justify the use of the non-relativistic approximation. 

2 Newtonian Approximation for Boson Stars 

We consider complex scalar fields which are coupled to gravity only. The 
action is given by: 

S = J d A x (-^ + g^dp&d,,® - M 2 $*$) (1) 

Since we will consider only the weak gravity limit of general relativity, the 
metric is expanded as g^ = rj^ + with \h^ u \ <C 1 and 
f]p, v = diag(l, — 1, — r 2 , — r 2 sin 2 (#)). Using the weak field approximation of 
General Relativity [8], we get 

□V = -16irGS^ (2) 
where = - \r]^T and 

Tpv = d^dv® + d^*d^ - rj^(rj a ^d a ^*d^ - M 2 $*$) (3) 
The equation for $, derived from (1), is: 

□$ + m 2 $ = (4) 

We will search for solutions with stationary rotation so $ will be allowed 
to depend on t and (p only through a phase, <£>(f, t) = (f)(r,8)e twt e imip . T^ U) 
given by (3), will be independent on both t and <p, and so will be the metric. 
To take into account deformations produced by the rotation, we allow both 
and h^y to depend on r and 6. 
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Finally, we will restrict ourselves to cases in which special relativistic 
effects are not important. The constraint imposed by this restriction will be 
checked to be consistent later on this paper. In this approximation, the only 
relevant component of is h 00 = 2V(r,9), where V(r,9) is the Newtonian 
potential. Equations (2) and (4) become : 

VV = MGM 2 (j) 2 (5) 

- D 2 (f) - (w 2 - M 2 )(f) + 2w 2 V(j) = (6) 

where D 2 (j) = (V 2 - V = V(r,6) and = <f>(r,9). In the non- 

relativistic limit, the gravitational binding energy E per particle must be 
much smaller than M, and the scalar field frequency may be written as 
w = E+M with \E\ M. The scalar field equation reduces to a Schrodinger 
equation: 

-^V 2 <P + MV<P = E<t> (7) 

Equations (5) and (7) must be solved subjected to the charge conservation 
constraint. The gauge invariance of the complex scalar field implies the 
conservation of j M = «(<9 M 00* — 0<9 M 0*) with conserved particle number: 

N = 2M J (f) 2 r 2 dr sin(0) d9 dip (8) 

3 Stationary Solutions 

We now look for solutions of (5), (7) with non-zero angular momentum de- 
scribed by $ = e lwt e imLp 4>(r,9) and V = V(r,9). We expand (f)(r,9) in asso- 
ciated Legendre functions: 
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1 oo 

and, since V, differently from $, has no ip dependence, we consider: 



(9) 



V(r,0) = Y,V l (r)P l (9) 



(10) 



With (9), (10) and the orthogonality relations of P™, we rewrite (5), (7) as a 
larger system of equations, which contains one equation for each value of I: 

K + ~K ~ ^f^o = GM\2k + 1) £ A lvl0 R lRl , (11) 

TV , j, 



2M 



2 jogo + 1) 



+ PP, 



■io 



oo oo 



(2 io + l)(i - ro ) 

2 (io + m) , ,, 



with ' = <9 r and 



(12) 



A lvl0 = £ iiffi WO P (*) (13) 
This system of equations may be rescaled by introducing the new variables 



[4]: 



f = NMr, 



V(r,0) = N 2 V(f,6), 



R(r) = N 2 (2G)- 1/2 R(f), 
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E = MN 2 E 



and 



N = GM 2 N 



(2/ + l)(/-m)! 



(l + m)\ 

With this new variables, our system can be written as: 



(14) 



' ' z i,i' 



(15) 



~ 2~, Z (Z + 1) - 

K lo + f K h f2 Kl ° 



+ ER, 



OO CO 



(2Z + l)(Z -m)! x . x . - f 

v u ' t=m I =m 



(16) 



where ' now stands for and we must add the normalization condition 

or 

derived from (8): 



/ 



Rf o (f) f 2 df = 1 (17) 

We may now systematically search for solutions of (15), (16) for different 
values of / and m. The basic idea is to consider solutions with R L ^ for 
one particular value of /, L, with Ri for I ^ L. The ground state equations 
correspond to the choice m = with Ri = for I ^ 0, in which case the 
equations simplify to: 

v 1 " + -v 1 ' = Ia 001 (r ) 2 



R'n + T-Rn 



+ ERq — - A 00 ;/_R V5 



2 ^ 

Z l'=0 



(18) 
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and A 00 i = j\ dx Pi(x) = 25i o. For I ^ 0, the equation for Vi is homoge- 
neous and the solution is the trivial Vj(r) = 0. In the expansion (10), the 
only non-zero component is V and we are left with the simple system: 

Vo + \% = (A,) 2 



2 - ' 



+ M = R Vo (19) 



The excited states will correspond to other choices of Rl ^ 0. As an 
example, we may focus on the I — 2, m — state described by the system 
of equations (15), (16) with the appropriate values of I and m. We set 
Ri = for I ^ 2. The r.h.s of equations (15) are given by ^^-A 2 2i ( y R2(f)) 2 
where ^22^ = /-1 dx (P2(x)) 2 Pi (x) , with the following numerical values: 
A 220 = 2/5, A 22 2 = 4/35, A 22 4 = 4/35, A 22 i = for / > 4 and for odd 
values of /. So, for I — 2 and m — 0, the Newtonian potential will be given 
by: 

V(r, 6) = V (r) + V 2 (r) P 2 {B) + V 4 (r) P 4 (6) (20) 

For / > 4, Vi, as solution of an homogeneuos system, may be set equal to 
zero. Vo, V2, V4 and i?2 will be given as solutions of: 

V " + -V^ = Ur 2 ) 2 (21) 
r 5 

V^ + ^-^ 2 = H(^ (22) 

v: + lv;- 2 p 4 = f 5 (R 2 y (23) 
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2 ~ 6 - 

R 2 + ~R 2 ~ ~^R 2 



+ ER 2 = R 2 



V + 2 -% + 2 -% 



(24) 



with the appropriate boundary conditions. For V , V 2 and V4, we may incor- 
porate these boundary conditions by formally solving (21), (22), (23) with the 
help of Green's functions. As usual, we consider the solution of: 



G' l >(r,r') + -G' l (r,r') 



-8(r — r') 



(25) 



given by : 



G,(r,r') 



(26) 



(21 + 1) r l + l 

which is regular at r — > and goes to zero for r — > 00. 

In the particular case we are considering, namely, / = 2 and m = 0, we 
end up with: 



Vc(r) = -- 



1 /- r°° 

- dr' {r'f {R 2 {r')f + dr' r'{R 2 {r')f 

T JO Jf 



(27) 



W = -35 



1 rf A /-co 1 

- / dr' (r') 4 (i? 2 (r')) 2 + r 2 / dr' - (i? 2 (r')) 2 
r io r 



(28) 



Mr) = ~ 



35 



1 rf „ /■oo I 

- / dr' (r'f (i? 2 (r')) 2 + f 4 / dr' _ (i? 2 (r')) 
r° jo Jr (r'J' 3 



(29) 



and we are left with just one (integral-differential) eigenvalue equation for 
R 2 (f): 



R 2 + —R' 2 — -^R-2 



R 2 W{f) = -E 2 R 2 



(30) 
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with 



W(f) = V (f) + 2 -V 2 {f) + ^V 4 (f) , (31) 
/ °° dff 2 {R 2 (r)) 2 = 1, and R 2 (0) = 0. 

The procedure described in the I — 2, m — example may be applied to 

all excited states. We should note here that the infinite sums in the r.h.s. of 

r i 2 

(15), (16) will always reduce to finite sums because P™(x) is a polynomial 
of order 1%, and may be expanded as a sum of Legendre polynomials with 
I < 2l . So, using the orthogonality of the Legendre polynomials, for any 
fixed value of l , A iolol = for I > 2l . 

Without any intention to make a complete analysis of these excited states, 
we present the numerical results for I — 0, m — 0, the groundstate, and for 
the excited states 1 = 1, m = 0;/ = l, m = l;/ = 2, m — 0. Basically, 
we start with a trial function, R l=0 (f), satisfying the appropriate boundary 
conditions. The index % gives the iteration order. We construct W(f) using 
(31), or the equivalent expression for other states, and we use it to obtain 
a new R t=1 (f), with the corresponding eigenvalue E l=1 . The process is then 
repeated with R t=l , and so on, until convergence is achieved in the solutions 
and, up to a given precision, no effect is introduced by new iterations. Our 
results are plotted in the figures 1, 2, 3. Note that the excited states are 
qualitatively different from the ground state and this is why they can not be 
obtained as simple pertubations of the ground state. 

4 Conclusion 

Following the prescription outlined here, in principle, we may construct all 
the excited levels of the rotating boson star, and based on these states, pre- 
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dictions can be made about the spectrum produced by boson stars decaying 
from excited states to states with lower energy. However, we still must check 
the consistence of the non-relativistic approximation made at the beginning. 

Basically, we must require that the gravitational field is weak and that 
no relativistic speeds are present. The weakness of the field is verified by 
requiring that the mass of the star in the Newtonian approximation is only 
slightly affected by rotation. The total mass of each state, M T , is given by 
[4]: 

M 2 

M T = MN + NE = N{ Pl ™ ck )(l + EN 2 ) (32) 
The rotation of the star changes its mass by the relative amount: 

= M^ ted - M° = (E exated N 2 - E Nq) 

M° 1 + E N 2 1 ' 

where we are comparing states with the same total number N of particles, 
and N is related to N and N by (14). To be consistent with the Newtonian 
approximation, we should only apply the above results for boson stars with 
A < 1. Since \E \ ~ 1CT 1 and \E excited \ < \E \, A <C 1 is always satisfied 
by Nq <C 1. For Nq ~ 1CT 2 , the boson star has a mass of order 10 10 kg, too 
small to resemble a conventional stellar object. However, this number is not 
much different from the maximum mass of boson stars of free scalars in the 
ground state, which is roughly 10 11 kg. 

The restriction coming from the small velocities limit is more severe. To 
check for relativistic speeds, we assume that a global effective angular velocity 
VL e ff may be associated to the stationary excited state, with: 

L e / f = IQ eff and 1 = J P z 2 dV (34) 
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and p is the energy density given by T 00 . Taking into account that I > Iq 
where J is the moment of inertia of the ground state configuration, we have: 



L e / f > J n eff = \ < rl > Mq n eff (35) 

with < rl >= = J l T J Po r 2 dV. 

This effective angular momentum is then identified to the angular mo- 
mentum of the configuration obtained by direct integration of the angular 
momentum tensor: 

f N 
L z = J (xT 0y - yT ox ) = mw— (36) 

Taking m = I and comparing (35) and (36), we have: 

S fff < « » (37) 
and v max ~ 2 < r > fi e// satisfies : 

V < 3 Tr^ V G N ° (2/ + l)(/-m)! (38) 
A sufficient condition to guarantee the use of non-relativistic approximation 
is: 

<f> w ~ Z(Z + m)! 
3 77f7 ^ *° (2/ + l)(/-m)! <<C L (39) 
Now, w ~ wo, but < r > grows fast with the excitation levels and even faster 
grows the factor tJt^tt^ m - 

° (2t+l)(t— my. 

Our conclusion is that the non-relativistic rotation is possible only for 
very small bosonic objects (small values for N ) and, even starting with 
small N , only the first excited states will be well describe by the Newtonian 
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theory. For higher excited states and objects with iVo ~ 1, relativistic effects 
will be important and we can rely only on a fully relativistic calculation. 
These restrictive limits tell us that the numerical results of our calculation 
can not be applied to large, astrophysically important stellar objects. Even 
so, we showed that small compact bosonic objects rotate producing a discrete 
set of excited states and we set limits on the applicability of the Newtonian 
approximation. Since the introduction of self-coupling is known to increase 
the maximum mass, it would be interesting to study how the self-coupling 
changes the rotational states, a problem we hope to address in the future. 
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Figure Captions 



Fig. la. Ground state (I — 0, m — 0) function, R(f), normalized to 
J* °° dff 2 (R(f)) 2 = 1, for the last twelve iterations. Ei =0ra=Q = —0.16. 

Fig. lb. Newtonian potential for each iteration in the same case (/ = 0, 
m = 0). 

Fig. 2. Excited states also normalized to / °° dr r 2 (R(f)) 2 = I. Curve 
1: I = 1, m = 1, E 1A = -0.025 ; Curve 2: I = 1, m = 0, £i, = -0.0071; 
Curve 3: / = 2, m = 0, E 2i0 = -0.0013. 

Fig. 3. Newtonian potential for I = 2, m = 0, showing contributions 
from Vo(f), V2(r) and V<i(r). The potential well, with constant, non-zero 
value for r — > is Vo(r). The other two represent ^(r) and V4(r). 
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